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The close similarity between the hierarchies of multiple-point correlation functions for the 
diffusion-limited coalescence and annihilation processes has caused some recent confusion, raising 
doubts as to whether such hierarchies uniquely determine an infinite particle system. We elucidate 
the precise relations between the two processes, arriving at the conclusion that the hierarchy of 
correlation functions does provide a complete representation of a particle system on the line. We 
also introduce a new hierarchy of probability density functions, for finding particles at specified 
locations and none in between. This hierarchy is computable for coalescence, through the method 
of empty intervals, and is naturally suited for questions concerning the ordering of particles on the 
line. 

PACS numbers: 02.50.-r, 05.40.-a, 05.70.Ln, 82.20.-w 

Diffusion-limited coalescence, + A — ► A, and annihilation, A + A — » 0, on the line, have long been known to display 
anomalous kinetics (different from the mean-field reaction-limited regime) and to belong to the same universality 
class pj-[31j. In fact, the similarities run deeper than that, as the full hierarchy of multi ple-p oint correlation functions 
in the two processes, when expressed in different length scales, are identical 0, 0, |U |2l|, Yet there exist 
important differences between the two processes. Most conspicuously, the density function for the gap x between 
adjacent particles (the so called inter-particle distribution function, or IPDF) falls off, for x —> oo, as e~ ax for 
coalescence, but only as e~° x for annihilation []J, |23|, |2fj, |3l| . This raises the question whether the infinite hierarchy 
of multiple-point correlation functions uniquely determines a system of particles on the line [2^ • 

In this communication, we answer this question on the affirmative: the hierarchy of multiple-point correlation func- 
tions does provide a unique representation of an infinite particle system. Indeed, the IPDF can be computed from 
the multiple-point correlation functions using an inclusion-exclusion formula, as shown in Eq. @- The correlation 
functions of the two processes are simply not the same, despite the similarity upon the rescaling of space. The precise 
relationship between coalescence and annihilation, and its consequences, is thoroughly discussed. An alternative hier- 
archy of probability density functions, for finding particles at locations X\,...,x n and none in between ^ is introduced, 
and we show how to compute it explicitly, for coalescence, using the method of empty intervals |26t I27L l28t l29t l30t l3^| . 
This new hierarchy is naturally suited for answering questions concerning the ordering of particles on the line. It is 
further used for obtaining an exact new expression, as well a series of systematic approximations for the IPDF of the 
annihilation process |23ll3l|. 

Equivalence between coalescence and annihilation 

The n-point correlation functions (the joint probability for finding particles at positions x\,...,x n at time t, 
simultaneously) of diffusion- limited annihilation and diffusion- limited coalescence in one dimension seem very similar. 
Indeed, for suitable initial conditions, one has the exact result h m, mill: 

pT\x x , ...,x n ;t) = ^(%*{xi, ...,*»; t). (1) 

The relation Q has the following simple interpretation: To obtain a configuration of the annihilation process, with the 
correct weight, select a configuration of the coalescence process but then retain only half of the particles, i.e., randomly, 
and independently, retain (or discard) each of the original particles with probability 1/2. 

This interpretation follows from the well-known observation that the two processes may be realized simultaneously. 
Starting from a random initial configuration of particles on the line, tag each of the particles with probability 1/2 and 
run a diffusion-limited coalescence process: When two particles meet they coalesce into a single untagged particle if 
the parents are alike (either normal or tagged), and into a tagged particle if the parents are different. Clearly, the set 
of all particles in the system (tagged and untagged) represents a configuration of the coalescence process A + A — > A, 
while the subset of tagged particles follows the annihilation process A + A — ► . At any time t a particle is untagged 
if and only if it has an even number of tagged ancestors at time t = 0, while the ancestors of two different particles 
form disjoint sets. It follows that at time t particles are tagged with probability 1/2, independently from one another. 
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For the suitable initial conditions for which Q holds, the density of particles in the anihilation process is half the 
density in the coalescence process. It is tempting to rescale space in the coalescence process by a factor of two, to 
impose the same density in both processes. Defining yi = 2xi and effecting the change of variables p^ oal (xi, . . .x n ;t) — > 
P™ al (2/i, ■■■y n ;t), one gets 

p™\x 1 ,...,x n ;t)=p™ al (y 1 ,...,y n ;t). (2) 

This relation has caused some confusion, leading one of us (DbA) to erroneously conclude that since the distribution 
of particles in coalescence and annihilation are different (as evidenced, for example, from their different inter-particle 
gap distribution functions) it must be the case that the infinite hierarchy of multiple point correlation functions does 
not uniquely determine an infinite set of points on the line |22j . Eq. J2J does not mean, however, that the functions 
/5* nnl and /5™ al are identical: they are applied to different arguments. Indeed, the stretching of space by a factor of 
two and eliminating half of the particles at random are not generically equivalent |24j . An easy way to see this is by 
considering the effect of the two operations on a lattice of equally spaced particles: stretching yields a new lattice with 
double the gap between particles, whereas the random elimination of half the particles leads to a disordered array. 



Inter-particle distribution functions (IPDF) 

The hierarchy of multiple-point correlation functions determines the system uniquely. Some quantities, however, 
such as the IPDF, are notoriously difficult to obtain, using the p n . This quantity is readily available for the coalescence 
process, through the method of empty intervals, but not for the annihilation process. 

To obtain the IPDF, one needs to compute the probability density P (xi,X2) for finding particles at X\ and x 2 
and no particles in between. The density probability function of the gap x between particles is then Pq(0,x)/p |25| . 
Pq{x\,x 2 ) is, in principle, available from the p n . Indeed, 

P {xi,x 2 ) = p 2 (xi,x 2 ) - / dzx p 3 (x 1 ,z 1 ,x 2 ) + dzi dz 2 Pa(xi, z x , z 2 , x 2 ) (3) 

J x± J X\ J z\ 

On the right hand side of this equation, events with exactly one particle between x\ and x 2 are counted once by the 
first term but cancelled by the second term; events with two particles in between are counted once each by the first 
and third terms, but are cancelled by the second term, which counts these events twice (z\ representing the first or 
second intervening particle), etc. In this fashion, only the events with no particles between X\ and x 2 are accounted 
for at the end. 

More generally, one can write down a similar expression for Pk(x\, x 2 ) — the probability density for finding particles 
at x\ and x 2 and exactly k particles in between: 

Pk{xi,x 2 ) = £ (-l) B-fc (? Wsi.a*), (4) 



n>k 

where 



R n {xi,x 2 )= J "'J P2+n(x 1 ,z 1 ,...,z n ,x 2 ) dz 1 ■ ■ ■ dz n . (5) 

x 1 <z 1 <---<z rl <x 2 

In principle, all the p n can be computed explicitly, both in the annihilation and the coalescence process, but their 
actual expressions are complicated and © becomes impractical. Another possibility is to derive P anm from the P™ al 
using the correspondence between the configurations of both processes: To obtain a configuration of the anihilation 
process with particles at positions X\ and x 2 and nothing in between, we start with a configuration of the coalescence 
process with particles at X\ and x 2 and exactly k other particles in between. Retaining each of these particles with 
probability 1/2, there is a probability 1/4 that the particles at x\ and x 2 stay put, and a probability (l/2) fc to kill 
the k particles in between. Therefore, we have 



nanm 
M) 



(x u x 2 ) = ±Y,^ P r\xuX2). (6) 

fc>0 



This is but a special case of Eq. (23) in [3l]], already suggested by Derrida and Zeitak. 

Obtaining the P%° (xi,x 2 ) from the p™ al , remains still an impractical proposition. Fortunately, in the case of the 
coalescence process, the method of intervals allows for a more direct way. 
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New hierarchy for finding n sequential particles 

Dimension 1 is special in that one can meaningfully discuss the ordering of the particles on the line. Thus, instead 
of the traditional hierarchy of multiple-point correlation functions, we propose an hierarchy designed to keep track of 
the sequential order of the particles. 

Let w„(ii, . . . ,x n ; t) denote the probability for finding particles at x±, x 2 , . . . , x n at time t, but no other particles 
in the intervals (x\, X2), (x 2 , x 3 ), ■ ■ ■ , (x„_i, x n ). Clearly, the uj n determine a distribution uniquely. Indeed, 



p 2 {x x ,x 2 ) = ^2 11 UJ 2+n(xi,z 1 , . . . ,z n ,x 2 ) dz 1 ■ ■ -dz n , (7) 

n=Q x 1 <z 1 <---<z„<x 2 

P3(xi,x 2 ,x 3 ) = /"'/ u 3+k+ i(xi,yi,. ..,Vk, X2,zi,. ..,z h x 3 ) dyi---dy k dzi---dzi, (8) 

7.7 j j 



fe >' x 1 <y 1 < - <y k <x 2 <z 1 < - <zi<x 3 



and similar expressions for p m , m > 3. Thus, defining lui(x) = pi(x), the complete hierarchy of multiple-point 
correlation functions, {p n }, can be derived from the hierarchy of sequential particles, 

The io n are better suited to deal with questions regarding ordered sets of particles. A relevant example are the P k , 
which instead of the infinite sum in are now simply given by 



Pk(xi,x 2 )= J -J u 2+ k(xi,z 1 , . . . ,z k ,x 2 ) dzi ■ ■ -dz k , (9) 

x 1 <z 1 <---<z k <x 2 

and, in particular, Pq(xi,x 2 ) — lu 2 (xi,x 2 ). 

In the case of coalescence, the hierarchies {p n } and {uj n } may be derived explicitly through the method of empty 
intervals. Specifically, the method of intervals yields expressions for E n (x\ , ?/i , . . . , x n , y n ; t) — the probability that 
the intervals (x%, yi), . . . , (x n , y n ) be simultaneously empty at time t |30|. The two hierarchies are obtained as spatial 
derivatives of the E n , but evaluated at different points: 

gn 

v^l ? * * * j 1 ^) — ~q ~q -^n 5 Vl 1 • ■ • 1 %n ) Un \ ^) \yi—X2-,---,y n -i =x n ,y n =x n j (-^) 



dxi ■ ■ ■ dx n 



Pn\X\i • • • , X n ] t) — E n (xi , yi , . . . , X n7 y n , t) \y ± — x\ ,...,y n — x n ■ (H) 



Actually, oj n can be computed somewhat more cheaply, from E n _i rather than E n : 

gn 

W n (xi, ...,X n ]t) = ~ q Xi ._q x ^g -E n -l(xi,yi, ■ ■ - ,a;n-l)J/n-i;*)li/i=a!a.".,*»-i=x»- ( 12 ) 

We have written Eq. pO(l merely to showcase the beautiful symmetry between p n and u> n . 

The first few w£ oal computed for coalescence, in the long time asymptotic limit, using the empty intervals derived 
in EH!!, are ESI 



u^\ Xl ,x 2 ;t) = V^ 2 a 2 e-«-, (13) 
^ oa \x 1 ,x 2 ,x 3 ;t) = ^ p ^ 13 ( e -&-& -e~&), (14) 
tol oa \x 1 ,x 2 ,x 3 ,x 4: ;t) = v^7> 4 {£u(e~^ 4 - e~^ 4 ~ 2 ^ 3 + e'^-ZL-tL _ e -fi2~tL + e ~&-&-tL _ e -£i3-&) 

+V^(Ci463e" ? " _?223 - fra^e" 5 "^ 4 + £i2&4e _ &~&)erfc(f2 3 )}, (15) 



where p = l/y/2nDt is the long time asymptotic density of particles in the coalescence process, £y = (xj — Xi)/y8Dt, 
and erfc(a;) = (2/y / 7r) f°° exp(— it 2 ) dit is the complementary error function. 

Returning to the question of the IPDF in the annihilation process, all the w™ al can be computed, in principle, but 
obtaining the exact p anm from JHJ and JHJ is not an easy task. Nevertheless, this approach can be used to generate 
efficient approximations of the IPDF. For instance, for small inter-particle gaps x = x 2 — X\, one need only keep 
the first few terms in JBJ, since the probability of finding several particles in the gap becomes negligibly smaller as 
their numbers increase. Indeed, using just the first term in jnj) yields an expression that matches the exact result 
(Eq (43) in [3]]]) to order a; 3 ; the first two terms improve the match up to order x 7 , and three terms up to x 12 . 
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However, truncating © in this fashion, at any order, gives terrible results for large x, as it predicts a Gaussian decay, 
exp(— 7ra 2 ), instead of the correct exponential decay. 

For large values of x, the simplest assumption that the gaps between particles are independent, 

UJ 2+k\ x i^ z i^ z 2,- ■ ■ ,Zk,x 2 ) « z ■ (lo) 

p K 

leads to to an exponential decay of the IPDF, PQ nm (0, x) ~ 1.6777 exp(—1.2685x), at large a; (where distance is scaled 
so that the density of particles is equal to 1). The same result was already derived in [23j, using a similar assumption 
of uncorrelated particles, and compares favorably with the exact result of PQ nm (0, x) ~ 1.8167 exp(— 1.3062x) [3l| . 
Systematic improvements are achieved by taking into account more correlations. For instance, assuming 

coa l, v ^'(n, Zl , 22)^'^, Z 3 , ■ ■ - < a '(z^l, 3: 2 ) , , 

UJ 2+l\ X l> Z li Z 2,---i z l,X2) i (J-') 

P 

(for I odd; for I even, the last term in the product is uj c 2 oal (z u x 2 j), leads to P o anni (0,a;) ~ 1.728976 exp(-1.285339ar). 



Discussion 



In summary, we have elucidated the exact relation between diffusion- limited one-species coalescence and annihilation 
in one dimension. The precise meaning of the similarity between the respective hierarchies of multiple-point correlation 
functions, Eq. has been clarified: configurations of the annihilation process are obtained by random elimination 
of half the particles in the coalescence process. This, however, is not equivalent to the stretching of space by a factor 
of two. 

We have also introduced a different hierarchy of probability density functions: the Lo n (xi, . . . ,x n ;t) — the probability 
for finding particles at x\,...,x n and no particles in between, at time t. This new hierarchy capitalizes on the 
topological constraints special to one dimension, and is better suited for answering questions concerning specific 
numbers of particles. 

Both the traditional multiple-point correlation functions, the p n , and the u n can be derived from the distribution of 
empty intervals on the line, in a way that highlights their relationship, Eqs. 1)11) [I. 1)110. The two hierarchies determine 
an infinite system of particles on the line completely, and in particular it is possible to express the p n using the u> n , 
and vice-versa, albeit at the cost of resorting to infinite sums of unwieldy integrals. In this sense, the best situation 
occurs when the distribution of empty intervals is available, for it yields the p n and uj n directly. This is the case for 
the coalescence process, but not for annihilation. For annihilation the p n are known, but the oj n are harder to obtain. 
Derrida and Zeitak have obtained the IPDF (directly related to UI2) exactly [3l|. Here we presented an alternative 
approach, based on the relation between the oj n and p n , that provides useful approximations. 

The situation is not yet completely clear even for coalescence. For example, consider the probability of finding two 
particles separated by a distance x and having exactly k particles in between, P™ al (0,x). For small gaps, x <C 1/p, 
we obtain P™ al (0, x) ~ x a ^ k \ with a(k) = 1, 4, 8, 13, for k = 0, 1, 2, 3, respectively. We have not yet found a 
satisfactory way to predict a(k) for arbitrary k. It remains the subject of future studies. 
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